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THE HUMAN PUPIL SERVOMECHANISM 



INTRODUCTION 

This Study, performed on a desk-top- sized PACE® 
TR-10 general purpose analog computer, describes 
the simulation of pupillary response in the human 
eye to an increase in light intensity. When the eye 
is exposed to such an increase in light intensity, 
the iris contracts, thereby reducing the pupil area 
and regulating the amount of light falling on the 
retina. 



The prediction by analytical means of K c -- that value 
of K which causes an oscillatory condition--is 
complicated greatly by the presence of the ex- 
ponential factor, e~0« 188 While modified versions 
of common methods (such as the Nyquist diagram) 
can be used, an analog simulation of the system 
permits the investigator to determine K c quickly and 
directly. 

PHYSICAL SYSTEM 



SYSTEM ANALYSIS 

This phenomenon has been investigated by Stark (1) 
who found that the pupil responds in a frequency 
range from zero to about 4 cps in a manner which 
can be described by a transfer function such as 
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The physical system can be described by the sim- 
plified block diagram of Figure 1. 
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Figure 1: Simplified Block Diagram of Pupil 
Servomechani sm 



which relates the "pupil-area signal" to retinal 
flux. The quantity 0.16 is the zero-frequency mag- 
nitude; the exponential term represents a pure 
time delay of 0.18 sec. 

This system as described is stable, but if the 
zero-frequency magnitude is replaced by an ad- 
justable constant, K, there will be a value of this 
constant which will place the system on the verge 
of instability. The system then will be self-excited 
and will produce oscillations termed "induced 
pupillary hippus". 



SCALING AND MECHANIZATION 

Figure 2 shows the computer diagram, time scaled 
by a factor of 10 (i.e., the simulated system is 
ten times slower than the physical system). The 
time delay, e~^« 18 ^, is simulated by means of the 
well-known Pade approximation (2) 
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and also has been time scaled. 
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Figure 2: Scaled Computer Diagram 



The potentiometer and amplifier assignment sheets 
are shown in Figures 3 and 4. 
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Figure 3: TR-10 Potentiometer Assignment Sheet 
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Figure 4: TR-10 Anplifier Assignment Sheet 
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RESULTS 



By applying a step function corresponding to a rapid 
change in light intensity input and assigning a dif- 
ferent value of K for each run, one can quickly find 
the value of K c . In this problem, the value found 
was K c = 1.97. A pencil and paper analysis based 
on approximating e TS by the first three terms of 
its series expansion yields the result K c = 1.82, 
with an oscillation frequency of about 66 cycles per 
minute. 

EQUIPMENT COMPLEMENT 

The major pieces of equipment necessary for this 
simulation include: 10 Operational Amplifiers (7 of 
which are used as integrators) and 10 potentiom- 
eters. 

Notes: 1. Gains of 2, 12, 20, etc. can be obtained 
with parallel combinations of standard gains 
of 1 and 10. For example, a gain of 12 is 
achieved with two gains of 1 and one gain 
of 10 in parallel. Thus 
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2. To obtain K > 1, change the gain from 
1 to 10 of the channel of amplifier #7 
which is fed by potentiometer #8. 



3. This study was derived from a problem 
in "Principles of Electronic Instrumenta- 
tion by Lynch and Truxal, McGraw-Hill 
Book Company, Inc., New York, 1962, pp. 
682-685. 
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